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During  the  past  year  of  research  under  Air  Force  grant  AFOSR-80-0020  we 
have  finished  our  investigation  of  large-scale  turbulence  at  periods 
characterizing  the  meteorological  mesosoale,  and  we  have  now  completed  the 
first  part  of  a  study  of  the  validity  of  various  objective  analysis  sohemes. 
Both  of  these  projects  involve  the  analysis  of  data  from  the  NOAA  M-S-T 
(Mesosphere/  Stratosphere/Troposphere)  radar  located  at  Poker  Flat,  Alaska* 
The  details  of  the  system  have  been  desoribed  by  Gage  and  Balsley  [1978], 
Balsley  et  al.  [I960],  and  Balsley  and  Gage  [1980].  The  data  used  in  our 
projeot  consists  of  high  time  resolution  data,  one  profile  every  4  min.,  for 
six  different  heights  in  the  troposphere  and  lower  stratosphere,  averaged  over 
one-hour  intervals.  This  type  of  data  was  not  available  from  any  other  source 
prior  to  the  M-S-T  radar  and  has  provided  the  opportunity  to  oarry  out  some 
unique  investigations  of  importance  to  synoptio  meteorology. 

Investigation  of  Turbulence  at  the  Meteorological  Masoscale 

Forty  days  of  horizontal  wind  data  from  the  radar  were  used  to  investigate 
turbulence  in  the  upper  troposphere  and  lower  stratosphere  at  periods  between  2 
hours  and  40  days.  The  power  spectra  of  the  winds  showed  a  k**5/3  inertial 
subrange  at  time  scales  characteristic  of  the  mesosoale.  Various  arguments  oan 
be  used  to  show  that  for  the  time  soales  in  question  the  flow  is 
two-dimensional  (see  Appendix  I).  The  implication  is  that  a  source  of  energy 
exists  at  spatial  soales  less  than  100  km,  propagating  energy  up  the  spectrum 
toward  lower  wavenumbers.  A  candidate  for  this  souroe  of  energy  is  the  release 
of  latent  he*.t  associated  with  oonvection. 

The  results  of  this  part  of  the  study  whioh  has  now  been  completed  are 
described  in  detail  in  the  paper  reproduced  in  Appendix  I.  The  paper  has  been 
submitted  to  J,.  Atmos.  Sol,  and  has  gone  through  the  refereeing  process.  The 
paper  was  deemed  to  qualify  for  "probable  acceptance  after  revisions." 


Sinoe  the  observations  of  th®  meteorological  data  network  ar®  irrogularly 
spaced  and  of  variable  density  in  various  parts  of  th®  globe,  some  t®ohnlque 
has  to  b®  used  to  interpolate  the  observations  to  a  regular  grid  that  oan  be 
used  as  input  for  a  numerioal  foreoast  model.  Beoause  the  Interpolation 
problem  is  suoh  a  basio  part  of  numerioal  veather  forecasting,  sohemes  to 
accomplish  this  have  existed  as  long  as  the  numerical  models.  The  first 
objective  analysis  algorithm  to  gain  wide  aooeptanoe  was  that  formulated  by 
Cressman  [19591.  It  was  computationally  efficient  and  seemed  to  provide  good 
agreement  with  the  results  of  subjective  analyses  of  the  same  data  fields  and 
was  used  by  the  National  Meteorological  Center  for  many  years.  Sinoe  then  the 
Cressman  technique  has  been  expanded  upon  by  various  people  who  have  tried  to 
include  more  of  the  aotual  physics  of  the  atmosphere  in  their  particular 
objective  analysis  method.  Some  of  these  have  been  reviewed  by  Kruger  [1969], 
and  others  were  described  by  Otto-Bliesner  et  al.  [1977]. 

The  goal  of  all  objective  analysis  procedures  is  to  generate  a  value  of 
the  meteorological  variable  in  question  that  is  representative  of  the  true 
value  at  the  grid  point  baaed  on  observations  at  a  duster  of  points  around  the 
grid  point.  The  interpolated  values  are  also  required  to  be  in  a  balanced 
state  suoh  as  that  represented  by  the  geostrophic  wind.  The  seoond  condition 
is  imposed  to  avoid  generating  spurious  waves  once  the  model  is  initialized. 
Since  the  balanoe  condition  is  an  integral  part  of  objective  analysis,  it 
provides  a  oheok  on  the  validity  of  the  various  sohemes,  particularly  if  a  good 
estimate  of  the  balanoed  wind  is  available v 

Verification  of  the  validity  of  the  various  objective  analysis  methods  has 
always  been  difficult.  As  an  example,  we  can  consider  the  method  used  by 
Kruger  [1969]  who  was  interested  in  evaluating  the  differences  between  the 
methods  of  Cressman  [19593  >  Gandin  [19631,  and  the  Central  Analysis  Office  of 
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Canada  [Kruger  and  Aaaalln,  19673*  Eaoh  of  thaaa  sohemes  uaaa  a  weighting 
funotion  whioh  only  dapanda  on  radial  dlstanoa  from  tha  grid  point.  An  obvloua 
oholoa  for  tha  waightlng  la  tha  autoocrralation  funotion  for  tha  mataorologioal 
variabla  in  quaatlon.  In  faot,  Craaaman'a  and  Gendin1 s  aohamaa  uaa 
approximations  for  tha  autocorrelation  funotion,  and  tha  Central  Analyala 
Offioa  method  aotually  oaloulataa  tha  autocorrelation  from  paat  data.  To  teat 
the  validity  of  tha  different  methods,  Kruger  [1969]  compared  the  mapa 
generated  by  applying  the  aohamaa  to  mapa  generated  by  a  subjeotlve  analyala  of 
a  500  mb  height  field,  A  similar  technique  haa  been  uaed  to  teat  other  aohemea 
[see  e.g, ,  Otto-Blieaner  et  al.,  19773*  Of  oourae,  thla  merely  teata  tha 
ability  of  the  objective  analyala  routine  to  produce  an  analyala  in  agreement 
with  that  produoed  aubjeotively  by  the  meteorologiat .  It  la  not  a  teat  of  the 
ability  of  the  routine  to  produce  valuea  corresponding  to  the  real  state  of  the 
atmosphere.  In  faot,  Otto-Blieaner  et  al.  [19773  found  that  the  various 
aohemea  tended  to  agree  more  with  eaoh  other  than  with  the  subjeotlve  analysis. 

Due  to  the  large  uncertainties  in  the  winds  determined  by  the  radiosonde 
data  network,  it  haa  been  difficult  to  teat  the  objective  analysis  methods  by 
means  other  than  those  described  above.  Since  the  balloons  obtain  a  snapshot 
profile  of  the  winds,  they  are  likely  to  be  contaminated  by  oaoillationa  of 
relatively  high  frequency  compared  to  the  synoptic  soale  motions  of  interest  in 
the  numerical  models.  The  Poker  Flat  M-S-T  (Mesosphere/ 
Stratosphere/Troposphere)  radar  has  provided  an  independent  set  of  wind  data 
that  oan  be  used  to  provide  a  test  of  the  objective  analysis  procedures 
described  above.  The  radar  obtains  a  wind  profile  once  every  4  minutes.  By 
averaging  these  over  longer  time  intervals,  the  higher  frequenoy  oaoillations 
can  be  eliminated  from  the  data,  and  a  good  estimate  of  the  quasi-steady  wind 
is  possible. 

In  our  study  to  date  we  have  investigated  three  different  analysis 
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methods.  A  forty  day  pariod  from  February  28  to  April  5*  1979  whan  the 
horizontal  winds  measured  by  the  radar  are  available  was  used  together  with  the 
standard  radiosonde  data  for  00Z  and  12Z  during  this  period  to  study  the 
agreeaent  between  the  radar  winds  and  the  geostrophio  wind  predicted  by  an 
objective  analysis  of  the  observed  height  field.  The  distribution  of  stations 
and  their  looation  relative  to  Poker  Flat  are  shown  in  Figure  1.  The  station 
furthest  from  Poker  Flat  is  Yakutat  looated  764  km  from  the  site  of  the  radar. 
Barter  Island  is  572  km  away,  MoOrath  476  km,  and  Anohorage  455  km.  The 
olosest  radiosonde  station  is  Fairbanks  whioh  is  36  km  SSE  of  Poker  Flat.  An 
analysis  whioh  uses  five  stations  is  usually  considered  good,  especially  if  the 
stations  are  within  a  radius  of  1500  to  2500  km  of  the  grid  point  of  interest 
[Kruger,  1969]. 

The  algorithms  for  the  various  sohemes  oan  all  be  formulated  in  the  same 
way  as 


T  ^  r  *  ^  e  ft 

f  «  R  ♦  r  cff?  ♦  l  CT(S?) 

8  8  i«l  1  1  j-1  J  J  8 


(1) 


where  fg 1  is  the  interpolated  value  at  the  grid  point  g,  Rg  is  an  initial  guess 
of  the  value  of  fg*,  f*°  is  the  observed  value  of  f  at  observation  point  i,  and 
Sj°  is  an  observed  variable  that  can  be  used  to  infer  the  value  of  f.  For 
instanoe,  Sj°  oould  be  the  wind  field  if  f  is  the  height  field.  The  C's  are 
the  weights  given  to  eaoh  of  the  observed  values.  Different  objective  analysis 
sohemes  ohange  the  way  that  the  different  weights  or  initial  guesses  are 
specified,  but  otherwise  (1)  is  perfectly  general. 

For  all  our  calculations  the  objective  analysis  scheme  was  applied  to  the 
observations  to  provide  values  of  the  height  field  at  five  grid  points.  One 
grid  point  coincides  with  the  looation  of  Poker  Flat.  The  other  four  are 
distributed  so  that  they  are  looated  one  each  at  a  distanoe  100  km  north, 
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south,  ssat,  and  wsst  of  Poksr  Fist,  The  initial  guasa  for  tha  halght  field  at 
aaoh  of  tha  grid  points  was  found  from  tha  average  valua  of  tha  halght  ovar  tha 
forty  day  pariod  at  aaoh  of  tha  observation  pointa.  Tha  thraa  olosaat  atationa 
to  any  ona  grid  point  vara  uaad  to  oaloulata  tha  horliontal  gradiant  of  tha 
halght  flald,  and  this  valua  of  tha  gradiant  was  than  uaad  to  oaloulata  an 
expeoted  valua  at  tha  grid  point.  This  prooadura  falls  into  tha  oatagory  of 
using  ollaatology  aa  an  initial  guasa.  A  variant  on  this  la  to  uaa  tha 
foraoaat  output  from  a  nuaarioal  modal  aa  tha  initial  guaas.  Hovavar,  this 
option  waa  not  availabla  to  us  ainoa  wa  did  not  hava  a  modal  for  that  purpoaa. 

Although  va  hava  oarriad  out  an  analysis  of  tha  halght  fialda,  vhioh  la 
standard  in  meteorology,  there  la  ona  alight  variation.  For  tha  forty  day 
period  that  va  are  disouaaing,  the  radar  was  getting  a  fairly  continuous  stream 
of  data  at  six  or  aaven  heights  in  the  upper  troposphere  and  lower 
atratoaphere.  Since  tha  range  gataa  of  tha  radar  are  fixed,  it  was  not 
possible  to  vary  tha  height  of  the  measurement  to  correspond  to  a  standard 
pressure  level.  To  oiroumvent  this  problem,  the  pressure  at  the  height  of  the 
radar  measurement  was  oaloulated  from  the  Fairbanks  radiosonde  data.  The 
halght  corresponding  to  this  pressure  was  then  found  for  eaob  of  the  other 
observing  stations.  This  should  not  make  any  differenoe  in  oarrylng  out  the 
1 2-hourly  objective  analysis  since  this  is  independent  of  previous  values. 
This  can  only  affeot  the  determination  of  the  oliaatologioal  averages  for  the 
period.  Tha  average  beoomes  tha  average  halght  at  aaoh  of  tha  observing  points 
corresponding  to  the  pressure  at  a  fixed  height  ovar  Poker  Flat.  In  any  oase, 
the  agreement  between  the  measured  winds  and  the  geostrophic  winds  calculated 
from  the  height  field  improved  substantially  when  the  Initial  guess  field  was 
inoluded  as  will  be  seen  later. 

The  first  method  that  wa  investigated  was  that  developed  by  Cressman 
[1959]  and  used  at  the  National  Meteorological  Center  for  a  number  of  years. 
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Referring  baok  to  (1),  it  oorreaponda  to  tha  oaaa  in  which  Cja«0  and  ia 
given  by 


i  cr2-°2) 

n  (R4+D4) 


(2) 


whara  R  ia  the  radlua  of  influence  and  D  ia  the  diatanoe  between  the  grid  point 
and  the  obaervation.  The  factor  n  ia  the  total  number  of  atationa  uaed  in  the 
analyaia  for  one  individual  grid  point.  The  weighting  ia  aero  if  D>R,  The 
weighting  funotion  ia  ahown  aa  the  curve  labelled  CR1  in  Figure  2.  The  aeoond 
method  inveatigated  ia  that  represented  by  the  weighting  funotion  labelled 
QANDIN  in  Figure  2.  The  third  ourve  labelled  CAO  ia  the  autocorrelation  ourve 
calculated  from  paat  data  by  the  Central  Analyaia  Offioe  of  Canada.  The  real 
autocorrelation  beoomea  negative  beyond  1500  km  whereaa  the  Creaaman  and  Qandin 
ourvea  never  become  negative.  The  Qandin  ourve  la  a  Qauaaian,  and  the  Creaaman 
ourve  ia  aero  by  definition  beyond  the  radlua  of  influenoe. 

The  effeot  of  not  inoluding  a  good  initial  gueaa  of  the  height  valuea  at 
the  grid  polnta  oan  be  aeen  by  oomparlng  Figurea  3a  and  3b.  Both  repreaent  the 
reault  of  a  Creaaman  analyaia  with  a  radlua  of  influenoe  of  1 ,000  km.  The 
helghta  are  indioated  in  the  lover  left-hand  corner  of  eaoh  frame.  The 
agreement  between  the  calculated  and  meaaured  winda  ia  indioated  by  aigma  which 
ia  the  atandard  error  given  by 


a 


i  ®  -  4 

KTT 


(3) 


Here  n  ia  the  number  of  pointa  in  the  time  aeriea,  ia  the  meaaured  value  of 
the  wind,  and  v^  ia  the  geoatrophio  wind  oaloulated  from  the  height  field. 
Approximately  635  of  the  meaaured  valuea  will  lie  within  ±o  of  the  eatimated 
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value.  Fop  the  oaae  with  no  initial  gueaa  large  errors  ooour  in  the  sonal 
ooaponent  in  particular.  The  error  la  reduoed  by  a  little  over  8  a/ a  at  the 
worat  height  whioh  la  10.39  km  when  the  initial  gueaa  la  inoluded,  There  is  no 
algnlfloant  difference  in  the  aerldional  ooaponent  when  the  initial  gueaa  ia 
oaltted.  It  ia  evident  that  the  initial  gueaa  Juat  eatabliahea  a  baae  level. 
The  objeotlve  analyala  aoheae  then  determined  the  variation  about  that  level  on 
a  given  day.  Thus  the  errors  are  larger  in  the  sonal  ooaponent  which  haa  a 
larger  mean  wind  than  the  aerldional  ooaponent.  The  error  ia  alao  larger  at 
10.39  ha  whioh  oorreaponda  to  the  peak  in  the  mean  wind  aa  a  funotion  of 
height. 

Certain  general  features  whioh  are  evident  in  Figure  3b  are  noteworthy 
beoauae  they  per slat  in  all  the  varioua  analyaes  that  are  preaented  here.  The 
smallest  standard  error  la  at  the  lowest  height  for  whioh  we  have  data,  namely 
3.79  ka.  The  largest  error  oooura  at  8.19  km  and  10.32  km.  These  helghta  are 
alao  in  the  range  where  the  largest  mean  winda  ooour.  Finally,  the  error 
begins  to  deoreaae  with  increasing  altitude  in  the  stratosphere. 

The  smallest  error  variance  was  found  to  ooour  for  a  radius  of  influenoe 
of  750  km,  although  the  difference  between  using  a  radius  of  influenoe  of  1,000 
km  and  750  km  was  lesa  than  0.5  m/s  at  all  heights.  However,  the  difference 
between  a  radiua  of  influenoe  of  2,000  km  and  one  of  750  km  was  to  inorease  the 
error  by  1  to  3  m/s.  This  is  shown  in  Figure  3o. 

The  Gandin  method  with  its  Gaussian  weighting  gave  the  beat  results  when  a 
radiua  of  influenoe,  i.e.,  a  Gaussian  half-width,  of  500  km  was  used.  The 
result  is  shown  in  Figure  3d.  The  errors  were  consistently  larger  than  the 
best  oase  for  the  Cressman  analysis.  However,  the  difference  between  the  two 
was  inalgnifioant,  being  no  greater  than  0,4  m/s. 

For  the  results  described  above,  the  analysis  method  has  been  varied  and  a 
12-hour  average  of  the  radar  winds  has  been  used  as  a  comparison.  The  effect 
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of  increasing  the  averaging  interval  la  to  daoraaaa  tha  appop.  Thla  oan  b« 
saan  In  Figupa  3a  which  ahowa  a  4 8-hour  average  of  tha  radar  winda  oompared  to 
tha  Craaaaan  analysis  with  a  radlua  of  influanoa  of  1,000  km.  Coapariaon  with 
tha  12-hour  average  rasulta  ahown  In  Flgura  3b  ahowa  that  tha  arror  haa  baan 
raduoad  by  aa  much  aa  4  n/a  at  8.19  km  In  tha  sonal  ooaponant.  Tha  arror  haa 
not  baan  aubatantlally  raduoad  In  tha  narldlonal  ooaponant.  Aa  tha  averaging 
interval  inoreases  bayond  12  houra,  tha  arror  In  tha  narldlonal  ooaponant  ataya 
relatively  oonatant.  Tha  arror  In  tha  aonal  ooaponant  oontlnuaa  to  daoraaaa* 
Thla  oould  ba  explained  by  a  aarlaa  of  wave  perturbationa  propagating  on  a 
aonal  ourrant. 

The  actual  aohana  daaorlbad  by  Craaaaan  [1959]  Involved  ualng  the 
weighting  function  that  we  have  daaorlbad  above  together  with  a  aarlaa  of 
auooaaalva  a oar a  to  Improve  tha  eatiaate  of  tha  interpolated  grid-point  value. 
In  a  aarlaa  of  five  aucoeaalve  aoana  tha  radlua  of  influanoa  la  raduoad  by  20 % 
of  ita  original  value.  Thua  only  the  oloaaat  atatlona  influenoe  the 
oorrectlona  applied  on  the  laat  aoan.  Tha  raault  of  applioation  of  thla  method 
la  ahown  in  Figure  3f.  Although  the  method  produoea  quite  a  few  more 
oaoillationa  in  the  praaaura  field,  making  it  appear  aa  variable  aa  the  wind 
field,  the  overall  effeot  la  to  inoreaae  the  error. 

Aa  a  atandard  for  oomparlaon,  Figure  3g  ahowa  the  reault  of  oaloulating 
the  geoatrophlo  winda  ualng  the  obaarvad  heighta  at  Anohoraga,  MoGrath,  and 
Fairbanks.  These  three  stations  ware  ohosen  beoausa  they  are  nearest  to  Poker 
Flat.  The  winds  calculated  in  this  way  have  a  strong  oaoillation  between 
successive  12  hourly  values.  Therefore  a  3  point,  equally  weighted  moving 
average  waa  applied  to  the  radiosonde-derived  winds.  When  compared  to  the 
1 ,000  km  Cressman  analysis  shown  in  Figure  3b,  tha  error  at  stratospheric 
heights  is  aotualy  smaller  when  the  three-station  method  is  used.  Below  the 
tropopause  the  Cressman  method  gives  better  agreement. 


The  only  other  study  that  we  are  aware  of  that  has  investigated  the 
geostrophy  of  the  wind  is  the  one  by  Angell  [1959j»  He  used  data  from  a 
constant  level  balloon  (transosonde)  at  300  mb.  The  winds  were  derived  by 
using  the  change  in  position  of  the  balloon  over  8  hours.  The  curves  presented 
by  Angell  correspond  to  a  standard  error  of  approximately  12  m/s.  The  300  mb 
level  is  close  to  9-km  altitude.  At  8.19  km  and  10.39  km  the  standard  error 
for  the  Cressman  analysis  compared  to  a  12  hour  average  of  the  radar  winds  was 
10.69  and  10.32  m/s  respectively.  Although  these  values  are  smaller,  they 
compare  favorably  with  the  values  calculated  by  Angell. 

We  have  found  that  there  is  no  significant  difference  between  the  Cressman 
and  Gandin  objective  analysis  schemes.  Although  the  magnitude  of  the 
ageostrophic  wind  component  increases  when  the  radius  of  influence  is  not  at  an 
optimum  value  (750  km  for  the  Cressman  method),  by  far  the  largest  deorease  in 
the  ageostrophic  component  results  when  tho  radar  winds  are  averaged  over 
longer  time  intervals.  Successive  scans  do  not  appear  to  give  an  improvement 
over  a  single  scan  with  the  Cressman  method  when  an  optimum  radius  of  influence 
is  used  for  the  single  scan.  Finally,  a  simple  geostrophic  wind  calculation 
using  three  radiosonde  stations  within  500  km  of  the  radar  was  competitive  with 
the  Cressman  analysis  when  compared  to  a  1 2-hour  average  of  the  radar  winds  at 
stratospheric  heights.  In  the  troposphere  the  objective  analysis  methods 
provide  a  better  estimate. 

The  results  described  here  are  not  only  important  as  an  evaluation  of  the 
Cressman  and  Gandin  sohemes.  They  also  have  implications  for  the  more 
sophisticated  objective  analysis  methods  described  by  Schlatter  [1975]  and 
Sohlatter  et  al.  [1976]  for  example.  In  general,  the  improvements  on  the 
simple  schemes  have  involved  using  both  an  autocorrelation  for  the  variables 
being  interpolated  and  also  crossoorrelations  with  the  other  observed 
variables.  Thus  the  objective  analysis  method  of  Schlatter  [1975]  uses  the 
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observed  winds  and  the  geostrophio  wind  relation  to  provide  additional 
information  for  interpolating  tKe  height  field.  Therefore  an  analysis  of  the 
geostrophy  of  the  wind,  suoh  as  the  one  described  here,  will  also  be  important 
for  the  more  oomplex  objective  analysis  methods. 

We  plan  to  describe  the  results  presented  here  in  a  paper  to  be  submitted 
to  Monthly  Weather  Beview.  The  second  part  of  our  study  will  involve  a 
comparison  of  the  radar  winds  with  the  analysis  method  of  Schlatter  [1975]. 
Finally,  the  estimates  of  the  ageostrophio  wind  components  derived  from  these 
comparisons  with  the  objective  analysis  methods  will  be  used  to  study 
geostrophio  adjustment  (see  Blumen,  1972).  We  intend  to  use  the  radar  data  to 
determine  if  inertial  waves  generated  during  periods  when  large  ageostrophio 
winds  are  present  can  be  detected. 
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Figure  2.  Weighting  functions  for  three 
different  objective  analysis  schemes.  CR1 
is  the  Creasman  method,  GANDIN  is  the  Gandin 
method  which  uses  a  Gaussian  weighting,  and 
CAO  is  the  Central  Analysis  Office  of 
Canada  autocorrelation  curve. 


«  t 


I 


3.  79  KM  51 G*  3.  79 


50  60  70  60  90  100  50  60  70  80  90  100 


Figure  3a.  Cressoan  analysis  for  the  period  from  February  25 
to  April  3,  1979.  The  radius  of  Influence  used  was  1,000  km. 
The  heavy  line  shows  the  geostrophic  wind  calculated  from  the 
analysed  data.  The  light  line  shows  the  radar  data  averaged 
over  12  hours  around  the  time  indicated.  The  horizontal  axis 
is  Julian  days. 


Figure  3b.  Same  as  Figure  3a  but  with  an  improved  initial 
estimate  of  the  height  field.  The  average  heights  at  each 
of  the  observation  points  during  the  forty  day  period  were 
used  to  calculate  an  estimated  value  at  each  of  the  grid 
points. 
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Figure  3c.  Same  as  Figure  3b  but  with  a  radius  of 
influence  of  2,000  km. 
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Figure  3d.  Comparison  of  the  Gendin  objective  analysis 
method  with  the  12-hour  average  radar  winds.  A  Gaussian 
half-width  of  500  km  was  used  for  the  weighting  function. 
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Figure  3f.  Cressman  analysis  with  five  successive  scans 
compared  to  12-hour  average  of  the  radar  winds.  During 
the  first  scan  the  radius  of  influence  is  2,000  km.  On 
each  successive  scan  the  radius  is  decreased  by  400  km. 
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Figure  3g.  The  geostrophic  wind  determined  by  using  the 
observations  at  McGrath ,  Fairbanks,  and  Anchorage  compared 
to  the  12  hour  average  of  the  radnr  winds.  A  three-point 
moving  average  has  been  applied  to  the  geostrophic  wind 
calculated  from  the  radiosonde  data. 
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Abstract.  Zonal  and  meridional  wind  measurements  mad a  with  tha  Pokar  Flat  ' 
MST  radar  ovar  a  40  day  pariod  ara  usad  to  oaloulata  tha  fraquanoy  power 
spactra  at  halghta  batwaan  5.99  and  14.89  km.  Tha  winds  used  in  tha  analysis 
are  1-hour  averages  of  samples  taken  every  4  min.  We  find  that  the  spectra 
follow  an  f~5/3  power  law  in  the  range  of  periods  from  2  to  50  hours.  If  the 
Taylor  transformation  is  valid  in  this  frequency  range,  this  would  imply  a 
k-5/3  wave  number  spectrum,  corresponding  to  an  inertial  subrange  for 
two-dimensional  turbulence  at  the  meteorological  mesosoale.  The  implication  of 
this  result  is  that  a  source  of  energy  exists  at  the  high  wave  number  end  of 
the  spectrum  with  energy  flowing  towards  lower  wave  numbers. 

1 .  Introduction 

The  study  of  atmospheric  turbulence  spectra  is  important  for  two  reasons. 
First,  the  shape  of  the  spectrum  shows  the  position  of  the  sources  and  sinks  of 
energy  and  enstrophy  and  how  nonlinear  interactions  distribute  these  quantities 
in  a  statistical  sense.  Seoond,  the  limits  of  predictability  of  the 
atmospheric  state  depend  on  the  form  of  the  wavenumber  power  law  as  shown  by 
the  work  of  Lorenz  [1969]  and  Leith  and  Kraichnan  [1972].  A  large  number  of 
studies  have  been  made  with  the  aim  of  defining  the  frequency  or  wavenumber 
power  laws.  In  general  they  fall  into  one  of  two  categories,  those  dealing 
with  large  spatial  and  temporal  scales  utilizing  radiosondes  [Kao  and  Wendell, 
1970;  Kao  and  Lee,  1977;  Julian  et  al.,  1970]  or  superpressure  balloons 
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[Mantis,  1970,*  Wooldridge  and  Reiter,  1970;  Desbois,  1975];  and  those  dealing 
with  small  scales  utilizing  data  from  instrumented  aircraft  [Kao  and  Woods, 
1964],  meteorological  towers  or  tethered  balloons  [Chernikov  et  al.,  1969]. 
The  part  of  the  spectrum  where  little  information  has  been  available  is  in  the 
meteorological  mesosoale. 

The  mesosoale  is  characterized  by  periods  between  a  few  hours  and  a  few 
days  and  spatial  scales  of  1000  km  or  less.  Meteorological  towers  or  tethered 
balloons  are  capable  of  resolving  frequencies  in  this  range  but  they  cannot 
operate  outside  the  planetary  boundary  layer.  Since  radiosondes  are  launched 
twice  per  day,  the  minimum  frequency  that  can  be  resolved  has  a  period  of  1 
day.  This  is  at  the  outer  edge  of  the  mesosoale  range.  The  superpres3ure 
balloons  which  are  flown  in  the  southern  hemisphere  have  the  potential  for 
investigating  subsynoptic  scale  motions  [see  e.g.  Cadet,  1978]  but  until 
recently  the  temporal  resolution  has  been  on  the  order  of  a  day  due  to  problems 
of  data  storage  and  transmission. 

The  NOAA  MST  [Mesosphere/ Stratosphere/ Troposphere]  rada"  located  at  Poker 
Flat,  Alaska  was  put  into  routine  operation  in  the  first  months  of  1979.  The 
system  has  been  described  in  detail  by  Balsley  et  al.  [1980],  Gage  and  Balsley 
[1978],  and  Balsley  and  Gage  [1980].  The  radar  uses  the  Doppler*  shift  of 
signals  back-scattered  by  turbulent  fluctuations  in  the  refractive  index  to 
measure  line-of-si ght  velocities.  By  using  three  beams  pointed  in  different 
directions,  the  complete  profile  of  the  vector  winds  can  be  determined.  The 
time  resolution  is  only  limited  by  the  integration  time  required  to  obtain  a 
good  signal-to-noise  ratio.  However,  the  fact  that  the  Poker  Flat  radar  is 
designed  to  operate  unattended  means  that  there  is  a  practical  limit  imposed  on 
the  time  resolution  in  order  to  avoid  changing  data  tapes  too  often.  For  this 
reason  the  time  between  successive  profiles  is  four  minutes.  In  1 979  when  the 
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data  sot  described  here  was  obtained,  only  a  quarter  of  the  full  system  had 
been  completed.  Thus  the  altitude  range  in  which  useful  measurements  could  be 
obtained  was  limited  to  heights  between  5.99  and  14.79  1cm.  This  time 
resolution  and  height  range  make  the  Instrument  ideal  for  investigating  the 
power  spectra  of  motions  at  synoptic  and  mesosoales. 

2.  Description  of  the  .Data  Set  and  Spectral  Analysis 

During  the  time  from  February  23  to  April  5>  1979  the  Poker  Flat  radar  was 
operated  in  a  mode  with  two  beams  at  15°  off- vertical,  one  pointed  roughly 
north  and  the  other  roughly'  east.  Computing  limitations  excluded  the 
possibility  of  operating  a  third,  vertically-pointing  beam  during  that  time 
although  this  situation  has  been  rectified  since  then.  For  this  original  data 
set  it  was  assumed  that  the  entire  contribution  to  the  line-of-sight  velocities 
was  due  to  a  horizontal  wind  vector.  The  antenna  consisted  of  a  100  m2  phased 
dipole  array,  one-fourth  the  size  the  array  will  have  upon  completion.  Sixteen 
100  kW  transmitters  were  used,  giving  a  peak  transmitted  power  of  0.8  MW  per 
beam.  This  configuration  together  with  a  1 5  ys  pulsewidth  provided  one 
complete  profile  every  4  min  up  to  an  altitude  of  14.69  1cm  and  a  height 
resolution  of -2.2  km. 

Coherent  integrations  and  the  Fast  Fourier  Transforms  of  the  received 
signals  are  done  on  line.  The  spectra  are  then  processed  at  the  Aeronomy  Lab 
in  Boulder,  Colorado.  The  processing  scheme  is  described  in  Carter  et  al. 
[1980]  and  provides  line-of-sight  Doppler  shifts  at  each  height.  The 
assumption  that  the  entire  contribution  to  the  line-of-sight  velocities  is  due 
to  a  horizontal  wind  vector  is  an  increasingly  good  assumption  as  the  averaging 
interval  is  increased.  One  way  to  understand  this  is  to  consider  the 
polarization  relations  for  gravity  waves  in  an  isothermal  atmosphere  (given  by 
Beer  [1974]  for  instance).  The  requirement  for  the  validity  of  the  hydrostatic 


4 


approximation  turns  out  to  be  that  «2«Wfc2f  where  w  is  the  angular  frequency 
of  the  motion  and  is  the  Brunt- Vaisala  period.  For  periods  of  half  an  hour 
or  more  the  contribution  from  the  vertical  velocity  component  is  oertalnly 
negligible.  However  for  samples  taken  at  intervals  less  than  this,  the 
vertical  velocity  component  can  contribute  significantly.  To  avoid  problems  of 
this  type,  the  4  min  data  was  averaged  over  one-hour  intervals  centered  on  the 
full  hour.  The  data  thus  oonsists  of  42  days  of  one-hour  average  horizontal 
winds.  For  this  study  40  of  the  days  were  used,  yielding  a  time  series  of  960 
points. 

The  time  series  was  treated  in  several  way3  before  the  spectra  were 
calculated.  First,  missing  data  points  were  replaced  by  interpolating  between 
adjacent  points.  The  nurd-ar  of  missing  points  was  less  than  45,  out  of  the 
total  of  960,  for  the  heights  discussed  here.  Above  14.69  km  and  below  5.99  km 
data  was  also  available,  but  the  number  of  missing  points  was  approximately 
four  times  as  large.  After  the  missing  data  points  were  replaced  by  the 
interpolated  values  the  mean  was  calculated  and  subtracted.  Then  a  Hanning 
window  was  applied,  affecting  96  or  10$  of  the  points  at  each  end  of  the 
series.  The  Hanning  window  helps  to  eliminate  energy  spillage  from  the  low 
frequency  end  of  the  spectrum  into  the  high  frequency  end  as  a  result  of 
broadening  due  to  discrete  sampling. 

The  original  data  is  shown  in  Figure  1 .  The  five  curves  at  the  top  of  the 
figure  are  for  the  zonal  wind,  and  the  curves  at  the  bottom  are  for  the 
meridional  component.  The  data  that  is  shown  has  not  been  treated  in  any  way. 
Both  the  average  for  the  entire  time  series  at  each  height  and  the  variance  are 
given  below  the  corresponding  curve.  Each  time  series  was  screened  for 
occurrences  of  wind  variations  greater  than  7  m/s  between  adjacent  points.  In 
almost  all  cases  where  variations  larger  than  this  occurred  a  large  change  in 
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the  wind  oooponent  was  found  both  preceding  and  following  the  questionable  data 
point.  A  data  point  such  as  this  was  considered  to  be  erroneous  and  replaced 
by  the  Interpolated  value.  Removing  these  erroneous  data  points  decreased  the 
amount  of  variability  in  the  high  frequency  end  of  the  computed  power  spectra. 

The  Fourier  transform  was  computed  with  an  algorithm  described  by 
Singleton  [19671*  This  routine  has  the  advantage  that  the  number  of  points 
does  not  have  to  be  a  power  of  2.  After  the  spectrum  was  calculated  it  was 
smoothed  by  the  funotion 

P(i)  a  1/64  (P(i-3)  +  6?(i-2)  +  15PU-1) 

+  20P(i)  +  1 5P(i+1 )  ♦  6P(i+2)  +  P(i+3) )  U) 
as  suggested  by  Endlich  et  al.  [19691.  The  smoothing  decreases  the  resolution 
but  increases  the  statistical  stability. 

3.  FraweasY  apgqtra  and,  aiM  sariaacaa 

The  resulting  spectra  are  shown  in  Figures  2a  and  2U.  Figure  2a 
represents  the  spectra  for  the  zonal  component  at  the  five  heights,  and  Figure 
2b  shows  the  corresponding  spectra  for  the  meridional  component.  Spectra  at 
succesive  heights  have  been  multiplied  by  a  factor  of  100  to  separate  them  on 
the  graph.  The  units  for  the  two  spectra  at  5.99-tan  altitude  are  correct. 
Least-squares  fits  of  a  funotion  of  the  form 

P(f)  a  P0  (f/f0)n 

were  made  to  determine  the  spectral  index  n.  Here  ?  is  the  power,  and  f  is  the 
frequency.  The  calculated  value  of  n  is  shown  next  to  the  corresponding  curve 
in  Figures  2a  and  2b.  A  dashed  line  is  used  to  show  the  fit.  The 
least-squares  computation  was  carried  out  for  frequencies  in  the  range  from 
f  a  0.015  hr"1  to  f  *  0.45  hr~1 .  The  average  value  of  the  spectral  index  for 
the  zonal  component  is  -1  .604  ±  0.279.  The  average  value  of  the  index  for  the 
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meridional  component  is  *1.604  ±  0.259.  The  value  of  *1.504  is  very  close 
to  -5/3.  There  is  no  particular  reason  to  expect  an  f-5/3  power  law;  however, 
if  the  Taylor  transformation  is  valid  in  this  range  of  frequencies,  then  a 
k-5/3  power  law  is  implied.  This  would  correspond  to  an  inertial  subrange. 
There  is  no  independent  evidence  from  this  study  to  indicate  whether  the  Taylor 
transformation,  i.e.,  that  temporal  scales  are  related  to  spatial  scales  by  a 
constant  faotor  of  the  mean  wind  over  the  sampling  tir  is  valid,  however, 
the  results  of  a  study  by  Grown  and  Robinson  [1979]  in  which  European 
radiosonde  data  was  used  were  that  the  Taylor  transformation  can  be  used  for 
spatial  soales  smaller  than  1200  Ian  and  at  least  down  tc  200  km,  the  lower 
limit  of  the  resolution  of  the  network  they  used.  The  corresponding  range  o 
periods  where  the  transformation  was  valid  was  in  the  range  less  than  2  days. 

Gage  [1979]  summarized  the  diverse  evidence  for  the  existence  of  a  k“?/3 
law  inertial  range  at  the  meteorological  mesoscale.  He  also  used  the 
variability  calculated  from  a  7  hour,  high  resolution,  balloon  wind  measurement 
experiment  and  wind  measurements  made  with  the  Sunset  VHF  Doppler  radar  in 
Colorado  over  a  14  hour  period  to  show  that  the  variance  defined  by 

<JT  *  {[v(t)-v(t+t)]2}1/2  (2) 

follows  a  power  law  out  to  time  lags  of  at  least  3  to  5  hours.  This  is 

consistent  with  a  k“S/3  power  law  when  the  Taylor  transformation  is  valid  as 

shown  by  Gage  [1979]. 

Plots  of  the  variance  aT  are  presented  in  Figures  3a  and  3b  for  the  zonal 
and  meridional  wind  components  respectively.  Lines  with  slopes  of  +1/3  are 
superimposed  on  the  curves  for  reference.  Successive  ourves  have  been 
multiplied  by  a  faotor  of  ten  to  separate  them  on  the  graph.  The  scale  is 

correct  for  the  variance  curve  at  10.39-km  altitude.  Overall,  the  1/3  slope 
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gives  a  reasonable  fit  to  tha  curvas.  Howavar,  thara  is  a  pronounoad  "bump" 
with  a  peak  naar  25  hours  at  8.19,  10.39,  and  12.59  km  which  is  particularly 
notioeable  in  tha  zonal  wind  component.  This  peak  in  tha  variance  should  not 
be  confused  with  a  diurnal  effect  since  a  peak  in  the  varianoe  at  25  hours,  for 
instance,  will  oorraspond  to  a  fraquanoy  with  a  period  of  50  hours.  This  can 
be  seen  in  tha  spectra  shown  in  Figures  2a  and  2b.  A  unique  feature  of  tha 
curvas  for  tha  meridional  wind  is  tha  rapid  decrease  in  the  variance  for  lags 
greater  than  600  hours  or  25  days. 

5.  Dlasuaalflfl 

In  the  preceding  section  it  was  shown  that  the  frequency  spectra  for 
motions  with  periods  in  the  range  of  2  hours  out  to  periods  around  50  hours 
generally  follow  a  -5/3  power  law.  Kao  and  Wendell  [1970],  Kao  and  Lee  [1979], 
Wooldridge  and  Reiter  [1970]  and  others  have  used  temporal  and  spatial 
transforms  of  radiosonde  data  to  show  that  at  the  synoptic  scale,  frequency 
power  speotra  follow  a  -1  power  law  while  wavenumber  power  spectra  follow  a  -3 
power  law.  Periods  of  2  to  3  days  are  at  the  high  frequency  limit  of  these 
studies.  Indeed,  it  may  be  argued  that  there  is  some  hint  of  a  change  to  a 
slope  of  -1_  at  periods  near  50  hours  in  the  spectra  shown  in  Figures  2a  and  2b. 
The  -1  slope  frequency  range,  corresponding  to  the  -3  wavenumber  regime,  is 
consistent  with  a  region  on  tha  high  wavenumber  side  of  a  source  of  enstrophy, 
presumably  tha  baroolinlc  Instability,  with  enstrophy  cascading  toward  larger 
wavenumbers  [Kralchnan,  1967].  Gage  [1979]  pointed  out  that  at  the 
meteorological  mesosoale,  characteristic  scale  sizes  are  too  large  for 
three-dimensional  turbulence  to  exist.  If  a  -5/3  inertial  subrange  is  present 
and  if  it  is  associated  with  two-dimensional  turbulenoe,  then  it  will  be 
characterized  by  an  energy  source  at  higher  wavenumbers,  with  energy 
propagating  up  the  speotrum  toward  lower  wavenumbers,  as  shown  by  Kralchnan 


8 


[1967]*  Sinoe  a  three-dimensional  inertial  subrange  is  also  characterized  by  a 
-5/3  power  law  but  has  the  energy  source  at  lower  wavenuabers  with  energy 
propagating  toward  higher  wavenuabers,  it  is  important  to  know  if  the  flow  is 
really  two-  or  three-dimensional.  Gage  [1979]  argued  that  the  validity  of  the 
hydrostatio  approximation  at  time  scales  corresponding  to  the  frequency  range 
studied  here  would  be  an  argument  in  favor  of  the  two-dimensionality  of  the 
flow.  We  can  also  consider  the  work  of  Charney  [1971].  He  was  able  to  show 
that  the  constraints  imposed  on  two-dimensional  turbulence,  i.e.,  conservation 
of  both  energy  and  enstrophy,  also  appear  if  the  flow  is  quasi-geostrophic. 
Therefore,  it  is  not  an  inherent  lack  of  variability  of  the  flow  in  the 
vertical  direction,  but  rather  the  geostrophy  of  the  flow,  that  makes  the 
two-dimensional  turbulence  theory  applicable.  A  study  by  Larsen  and  Kelley 
[1980]  using  the  same  data  base  discussed  in  this  paper  indicates  that  the 
geostrophlc  approximation  is  reasonable  for  winds  averaged  over  periods  of  half 
a  day  or  less.  This  result  also  argues  in  favor  of  the  applicability  of  the 
two-dimensional  rather  than  the  three-dimensional  turbulence  theory  for  flow  at 
the  meteorological  mesoscale.  Since  the  theoretical  work  of  Kraichnan  [1967] 
has  also  shown  that  enstrophy  cannot  propagate  through  a  region  characterized 
by  a  -5/3  power  law,  and  energy  cannot  propagate  through  a  region  characterized 
by  -3  power  law,  this  implies  that  there  must  be  a  sink  of  energy  and  enstrophy 
at  the  break  in  the  speotrum.  In  our  case  the  break  appears  at  periods  near  50 
hours.  Using  the  mean  value  of  the  zonal  wind  over  the  40  day  period,  this 
frequency  converts  to  a  scale  size  of  roughly  1800  km.  In  his  study  Gage  found 
that  the  break  oocurs  at  scale  sizes  around  1000  km.  The  physical  picture 

outlined  here  is  summarized  in  Figure  3  which  shows  the  various  power  laws, 

y 

direction  of  flow  of  energy  and  enstrophy,  and  souroes  and  sinks  at  their 
locations  relative  to  changes  of  slope  in  the  speotrum. 
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6.  CQOQluaiQn 

The  Poker  Flat  MST  radar  has  provided  a  -high  time  resolution  data  series 
of  horizontal  winds  over  a  40  day  period.  Speotral  analysis  of  these  winds 
indicates  that  the  two-dimensional  inertial  energy  subrange  discussed  by  Qage 
[1979]  does  exist  at  the  aeteorologlcal  mesosoale.  This  result  together  with 
results  of  earlier  studies  of  power  spectra  for  synoptio  scale  motions  implies 
that  souroes  of  energy  and  enstrophy  are  present  at  both  large  and  small 
soales.  The  large  scale  source  is  presumably  due  to  instabilities  of 
planetary- scale  waves.  At  the  small-scale  end  of  the  spectrum  the  energy 
source  is  characterized  by  periods  less  than  2  hours.  Theoretical 
considerations  imply  that  a  sink  of  energy  and  enstrophy  must  exist  at  the 
break  between  the  synoptic  and  mesosoale.  This  sink  is  characterized  by 
periods  of  approximately  two  days  and  scale  sizes  of  1000  to  2000  km. 
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Figure  1 .  This  shows  ths  original  data  sat  of  960  hourly  valuta  at  aaoh  of  tha 
fiva  half  fats  usad  in  tha  study.  Tha  avaraga  valua  and  tha  varianoa  at 
aaoh  of  tha  haights  art  shown.  Tha  top  hulf  of  tha  figura  shows  tha  zonal 

wind  component,  and  tha  bottoa  half  shows  tha  aaridional  oomponant. 

Figura  2a.  Thasa  art  tha  powar  spaotra  for  tha  tint  sarias  of  tha  sonal 
oomponant  shown  in  tha  pravious  figura,  Tha  paramatar  n  at  tha  right  of 
aaoh  spaotrum  is  tha  least-squares  fit  valua  of  tha  spaotral  index.  Tha 
avaraga  valua  of  tha  spaotral  index  is  -1.604  ±  0.279.  Spaotra  at 
suooessive  heights  have  bean  multiplied  by  a  faotor  of  100  to  separata 
them  on  tha  graph.  Tha  axes  are  correct  for  tha  spaotrum  at  5.99  km. 

Figura  2b.  Same  as  figura  2a  but  for  tha  meridional  wind  oomponant.  Tha 

avaraga  valua  of  tha  spectral  index  is  -1.604  ±  0.259. 

Figura  3a*  This  shows  tha  variance  of  zonal  wind  as  a  function  of  tha  lag. 
Tha  varianoa  is  defined  in  tha  text.  Curves  at  successive  altitudes  have 
bean  multiplied  by  a  faotor  of  10  to  separata  them  on  the  graph.  Tha 
vertical  axis  labeling  is  correct  for  the  varianoa  at  10.39-km  altitude. 

Figura  3b.  Same  as  figura  3*  but  for  tha  aaridional  wind  component. 

Figura  4.  This  graph  illustrates  tha  wave  number  powar  laws  that  are  discussed 
in  tha  text.  Tha  direction  of  tha  flow  of  energy  and  enstrophy  are  shown 
together  with  tha  various  souroas  and  sinks. 
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